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A characteristic feature of small- x lepton-proton data from HERA is geometric scaling — the 
fact that over a wide range of x and Q 2 all data can be described by a single variable Q 2 /Q^ix), 
with all ^-dependence encoded in the so-called saturation momentum Q S a.t(x). Here, we observe 
that the same scaling ansatz accounts for nuclear photoabsorption cross sections and favors the 
nuclear dependence Q 2 &t a oc A a Q 2 at , a ~ 4/9. We then make the empirical finding that the 
same A-dependence accounts for the centrality evolution of the multiplicities measured in Au+Au 
collisions at RHIC. It also allows to parametrize the high-p t particle suppression in d+Au collisions 
at forward rapidities. If these geometric scaling properties have a common dynamical origin, then 
this A-dependence of Qs ati A should emerge as a consequence of the underlying dynamical model. 



1. All data for the photoabsorption cross section 
tr 7 p (x,Q 2 ) in lepton-proton scattering with x < 0.01 
have been found Q to lie on a single curve when plot- 
ted against the variable Q 2 /Q 2 at , with Q 2 at ~ x~ x and 
A ~ 0.3. To further explore this empirical property of 
geometric scaling, we study here how experimental data 
on lepton-nucleus collisions constrain the geometric in- 
formation entering Q 2 at . We also ask to what extent the 
geometric scaling ansatz can account for characteristic 
features of particle production in other nuclear collision 
systems. 

Geometric scaling is usually motivated in the QCD 
dipole model where the total 7* ft cross section reads 

<r^(x,Q 2 ) = J d 2 vj\z\^ L (Q 2 ,v,z)\ 2 a'^(v,x).{\) 

Here ^t,l ar e the perturbatively computed transverse 
and longitudinal wave functions for the splitting of 7* 
into a qq dipole of transverse size r with light-cone frac- 
tions z and (1 — z) carried by the quark and antiquark 
respectively. Both, for a proton [h — p] and for a nu- 
cleus [h = A], a'^ ip (r,x) can be written as an integral 
of the dipole scattering amplitude Nh over the impact 
parameter b, 

a^ ip (r,x) = 2 J d 2 hN h {r,x-h). (2) 

In this setting, geometric scaling corresponds to the con- 
dition Nh(r,x;b) = N(rQ sat (x,h)). This can be seen 
by rescaling the impact parameter in Q in terms of the 
radius Rh of the hadronic target, b = b/ ^J-nR\, 

4 h L (x,Q 2 ) = nR 2 J ' d 2 vJdz\^ L {Q 2 ,v,z)\ 2 

x 2 J d 2 bN h {rQ sat (x,b)). (3) 

Since L (Q 2 ,r, z)\ 2 is proportional to Q 2 times a 
function of r 2 Q 2 , Eq. © depends solely on r = 



Q 2 /Q 2 at (x, b). In the case of 7*-A interactions, geomet- 
ric scaling is the property that the A-dependence of the 
ratio <Ty ^/ttR^ can be absorbed in the A-dependence of 
Q sa t,A(a;,b), 

^' A (ta) o^ijp = t a ) m 
ttR\ ttR 2 1 ' 

For this A-dependence, we make the ansatz that the sat- 
uration scale in the nucleus grows with the quotient of 
the transverse parton densities to the power 1/6, 

QL,A = QL, P [^j ^r A =r h ^J , (5) 

where the nuclear radius is given by the usual 
parametrization R A = (1.12A 1 / 3 - 0.86A" 1 / 3 ) fm. We 
treat S and irR 2 as free parameters to be fixed by data. 
2. In Fig. ^ we plot the experimental j*p data Q 
with x < 0.01 as a function of r = Q 2 /Q 2 at - F° r 
Q 2 at , we use in this plot the Golec-Biernat and Wiisthoff 
(GBW) parametrization [5| with i? 2 , = l/Q 2 at = (x/x ) x 
in GeV" 2 , x Q = 3.04 • 10~ 4 and A = 0.288. To 
safely extend to low virtuality, the ^-dependence of the 
GBW parametrization is modified by a mass term x = 

x Q2 K/ ) with nif — 0.14 GeV. The data are seen 
to be parametrized well by the scaling curve 

a^P{x, Q 2 ) = $(r) = a [ lE + T (0, + lnf] , (6) 

where je is the Euler constant, T(0,£) the incomplete 
r function and £ = a/r h , with a = 1.868 and b = 0.746. 
The normalization is fixed by (To — 40.56 mb. The func- 
tional shape of 10 can be motivated by making simpli- 
fying assumptions in evaluating (Q. For our purpose, 
however, Eq. © is just a convenient ansatz for the scal- 
ing function <&(t). 

To determine Q 2 at A , we compare the functional shape 
of © to the available experimental data for 7* A colli- 
sions with x < 0.0175 jl 0, 0, using £ = a/r\. The 



2 



parameters S and irR 2 in Q - JBJ are fitted by x 2 min- 
imization adding the statistical and systematic errors in 
quadrature. The data sets @, and have addi- 
tional normalization errors of 0.4%, 0.2% and 0.15%; the 
quality of the fit improves by multiplying the data by 
the factors 1.004, 1.002 and 0.9985 respectively. We ob- 
tain S = 0.79 ± 0.02 and 7ri? 2 = 1.55 ± 0.02 fm 2 for a 
X 2 /dof — 0.95 - see Fig. \I\foi comparison. If the nor- 
malizations are all set to 1, we obtain an almost identical 
fit with S = 0.80 ± 0.02 and 7ri? 2 = 1.57 ± 0.02 fm 2 for a 
% 2 /dof = 1.02. If we impose 6 = 1 in the fit, which cor- 
responds to Q 2 at oc A 1 / 3 for large nuclei, a much worse 
value of x 2 /dof = 2.35 is obtained. We conclude that 
the small- a; experimental data on 7* A collisions favor an 
increase of Q 2 at A faster than A 1 / 3 . The numerical coin- 
cidence b ~ S is consistent with the absence of shadowing 



in nuclear parton distributions at Q 3> Q 
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FIG. 1: Geometric scaling for 7*p (upper panel, data 
from Q), 7*^4 (middle panel, data from [M E|) and the ra- 
tio of data for 7* A divided by the scaling curve © (lower 
panel). Also shown in the lower panel are the data from Q 
for ratios over C. 

3. Can geometric scaling, and in particular the A- 
dependence and energy dependence of Q sa ,t,A(x), account 
for the pt-integrated multiplicity in symmetric nucleus- 
nucleus collisions at mid-rapidity? To address this ques- 
tion, we turn now to the heuristic ansatz 



dN 



AA 



dy 



sat, A 



ttR, 



(7) 



which arises in several models of hadroproduction 0, 0, 
ll(i Hl| . These models relate the parton distribution 
measured in a 1 A to the hadroproduction measured in 
nucleus-nucleus collisions. For example, the factorizcd 
formula |8| calculates gluon production by convoluting 
yl-dependent gluon distribution functions 



dN AB as 
oc — 

dyd 2 p t d 2 b p 2 



d 2 k <j> A (y, k 2 , b) <p B (y, (k - p t ) 2 ,b) , 



where 4>h{y , k , b) = J d 2 r expjir • k} Nh(r, x; b)/(27rr 2 ) 
For geometric scaling, c/> A (y, k 2 ,b)= <j){k 2 /Q 2 at A (y,b)), 
we find the dependence of Eq. (JZJ) , 



dN AA 



dy 



oc 



?sat,A 7r ^i 



pi 

d 2 S 



dk 2 d 2 b 




d 2 rd 2 b 



(9) 



Also without invoking factorization in Q , an y in tegrand 
with (/c/<2sat,A)-scaling leads to Eq. J7J), see [lfllll|. In 
all these models, a one-to-one correspondence between 
parton and hadron yields is assumed. 

To write (Q) in measurable quantities, we express the 
energy dependence of the saturation scale in terms of 
the GBW parameter A = 0.288, and we translate its 
A-dependence to an Ap art dependence fixed by our fit 
parameter 5 = 0.79 ± 0.02, 



1 dN 



AA 



N 



part 



drj 



N y 



1 tfplS 



(10) 



The overall normalization, independent of the energy and 
the centrality of the collision, is fixed to No — 0.47. As 
seen in Fig. |2 this reproduces without further adjust- 
ment experimental data from the PHOBOS Collabora- 
tion ^3] on charged multiplicities in Au+Au collisions at 
Vs = 19.6, 130 and 200 GeV/A. Even the pp data (0, 
as quoted in 01) at V$ — 19-6 and 200 GeV are ac- 
counted for by Eq. l|10fl . In the same figure, we show the 
result of (fTTTf) for intermediate RHIC energy (y/s — 62.5), 
for LHC energy (5500 GeV/A) and for smaller colliding 
nuclei. Eq. (|10fl implies that the energy and the central- 
ity dependence of the multiplicity factorize, in agreement 
with the results by PHOBOS [3. 

4. In the current debate of RHIC data on the sup- 
pressed high-pt hadroproduction in nuclear collisions, the 
relevance of nuclear shadowing has been discussed re- 
peatedly [3 HU- It is clear by now [n| that the A- 
dependence of ^-differential hadroproduction in nucleus- 
nucleus collisions and in deuteron-nucleus collisions at 
mid-rapidity both involve additional nuclear effects which 
are at least as significant as nuclear shadowing. On the 
other hand, arguments have been put forward 0, 0] 
that in d+Au collisions at forward rapidity, nuclear shad- 
owing may be the dominant effect. Motivated by the 
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FIG. 2: Energy and centrality dependence of the multiplic- 
ity of charged particles in AuAu collisions (| 1 01 compared to 
PHOBOS data |T^| . Also shown in the lower panel are the pp 
data [H| and results for y/s = 62.5 and 5500 GeV/A. 



phenomenological success of the scaling ansatz (|1(J|) . we 
now test to what extent the centrality dependence of 
^-differential hadron spectra in d+Au emerges natu- 
rally from the geometric scaling found in 7* A We start 
from the model JSJ, using for the fc t -diffcrential nuclear 
gluon distribution the scaling function ^(t) at virtual- 
ity $(r) ~ ct>(k = Q/2), such that r = fc 2 /4Q 2 atA , 
OLt A = ^cQs a t a/Cf- This approximation can be mo- 
tivated in a momentum space representation H3 of p. 
The parton distribution in the deuteron is taken to fall 
off sufficiently quickly, ~ l/fc™, n ^> 1, so that we can 
write 



dN dAu 

N coUl dr]d z p t 
dW dAu 

52 

N ca u 2 dr)d 2 p t 



Ncolli^CPi/Qsata) ^VcoUi^(T2) 



(11) 



We see the use of this rough pocket formula mainly in 
emphasizing the plausible claim that the suppression of 
d+Au at forward rapidity traces directly the suppression 
of nuclear parton distributions at small x. For the com- 
parison m Fig. El to data on the normalized yields of 
central and semi-central over peripheral dAu collisions, 
we use the number of collisions N co \\ in different central- 
ity bins [3 13.6 ± 0.3, 7.9 ± 0.4 and 3.3 ± 0.4. Only 
the two most forward rapidities 77 = 2.2 and 3.2 are com- 
pared. We find that Eq. ifTT|) captures main features of 
the recent data by BRAHMS [18J but it shows a weaker 
rapidity dependence. A more quantitative discussion is 



certainly beyond the accuracy of The only con- 

clusion from this exercise is that the more differential 
analysis of JHJ) is not inconsistent with data in d+Au. 

5. We now comment on the differences with other ap- 
proaches. The geometric scaling in 7* A data has been 
studied in 0, where a growth of Q 2 at A ocA",a<l/3 
has been found. This disagreement with our finding could 
have several origins. First, 0.01 < x < 0.1 was allowed 
in [J^; however, in this antishadowing region, we find no 
scaling in the data, as expected. Moreover, [ijj does not 
modify the variable x for small-Q 2 as done in this work 
and in the GBW model. Second, uses Ra oc A 1 / 3 
which leads to differences in particular for small A - we 
find a much worse fit in terms of x 2 for such an ansatz. 
In ^1 this dissagreement is improved by introducing 
a free parameter 7 in the Adependent normalization of 
the nuclear data, A -7-1 . In our case, however, this 
normalization is fixed by a dimcnsionful quantity given 
by the scaling condition 

The multiplicities in Au+Au collisions at RHIC have 
been studied jlo| on the basis of Eq. © assuming 
OLt A 00 A 1 / 3 . These authors are lead to an expres- 
sion which is Eq. (|10|l with <5 = 1 times an additional 
factor In (y / s A -^part) argued to come from scaling vio- 
lations. We note that for the accessible range of A, 
A 4 / 9 ~ A 1 / 3 ln(A/ 3 ) - this is the reason why both ap- 
proaches provide a fair description of the data at RHIC. 
However, the energy dependence in the logarithmic pref- 
actor introduced in |20| implies a flatter centrality depen- 
dence with increasing energy - this difference to Eq. i|10|) 
becomes sizable at higher t/s. 

Finally, the connection between the small x- and A 
dependence of parton distribution functions, and the sup- 
pression of normalized yields in d+Au collisions at 
forward rapidity has been discussed in several recent 
works [3 [III . Eq. JTTJ contributes to this discussion 
by illustrating to what extent the suppression of high-p f 
particles in d+Au at RHIC can be accounted for by the 
shadowing in 7 *A collisions ($ = a~t* A ). 

6. Here, we have discussed to what extent data for differ- 
ent collision systems in 7* A, dA and AA can be related 
through geometric scaling. Our study does not exclude 
the possibility that geometric scaling in j*p and j*A 
is a numerical coincidence without any dynamical ori- 
gin. However, geometric scaling also emerges naturally 
in non-linear small- a; QCD evolution equations |2jJ, |22| , 
which allow to absorb the entire dependence of small- a; 
parton distributions on energy and geometry into a sin- 
gle dimensionfull quantity, Q sat . The data discussed here 
are currently considered |l(j to provide the main support 
for such non-linear saturation effects. In fact, the scal- 
ing function $ in © resembles the asymptotic solution 
of the Balitsky-Kovchegov (BK) equation: it behaves 
as ln(fc/Q sat ) [(QLA 2 )"] for small [large] k [3 E^. 
Given that these non-linear evolution equations hold in 
a novel high-density regime of QCD which may become 
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FIG. 3: Normalized ratios of central and semi-central to pe- 
ripheral dAu collisions measured by BRAHMS [lsf compared 
to results from Eq. 1111 . The bands represent the uncertainty 
in the determination of N co n [Tg|. Results for the same cen- 
trality classes at the LHC are given in the lower panel. 

experimentally accessible, it is of obvious interest to ask 
whether the connection between geometric scaling in the 
theory and in the data can be made more quantita- 
tive. On the theoretical side, this requires at least the 
study of the impact parameter dependence [24| of small- 
x evolution, and the control of higher order effects. In 
particular, running coupling effects are known qualita- 
tively to decrease the energy dependence [25| and the 
A-dependence j26| of the saturation scale in comparison 
to the BK equation at fixed coupling (e.g., Refs. [26| fa- 
vor a weak A-dependence Q^ at A cx A a , a <C 1/3 but do 
not address the impact parameter dependence which is 
expected to increase the A-dependence). While an A- 
dependence of Ql at a x j 4 1 ^ 3 is often assumed [2^. much 
stronger ones (such as a ~ 2/3 |28|) have also been pro- 
posed. The present work has analyzed to what extent 
data constrain these energy- and ^-dependences. These 
constraints have to be met by non-linear small- a; evolu- 
tion or by any other model which aims at providing the 
common dynamical origin for geometric scaling in differ- 
ent nuclear collisions. 
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